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Abstract. The enhancement factors of the electric dipole moment (EDM) of the ground states 
of two paramagnetic atoms; rubidium (Rb) and caesium (Cs) which are sensitive to the electron 
EDM are computed using the relativistic coupled-cluster theory and our results are compared 
with the available calculations and measurements. The possibility of improving the limit for 
the electron EDM using the results of our present work is pointed out. 



1. Introduction 

With excellent progress in the experimental techniques and the computational advancements in 
the field of high precision atomic physics in recent years, it has now become possible to undertake 
rigorous studies to test the validity of many models of particle physics which question the very 
fundamental laws of physics [1, 2]. One such challenge for both theorists and experimentalists 
is the simultaneous existence of parity (P) and time reversal (T) violations giving rise to the 
electric dipole moment (EDM) of an elementary particle like the electron [3, 4]. It is interesting 
for theorists because it is a direct proof of the time reversal violation in nature and also it has 
profound consequences ranging from sub-atomic physics to cosmology [5, 6]. It will be a break 
through for experimentalists as it requires considerable ingenuity and the state of the art tech- 
nology to measure the tiniest electromagnetic moment to date. 

An atom can have EDM even in the absence of any electromagnetic field due to an intrinsic 
EDM of any of its constituent particles like electrons, nucleons, quarks or also due to P & T 
violating electron-nucleus interactions. In this paper, we have restricted ourselves to the cal- 
culation of the atomic EDM arising from the electron EDM. The intrinsic EDM of an electron 
interacts with the internal electrostatic fields in an atom, thereby enhancing the atomic EDM. 
Following the ideas first discussed by Schiff [7], Sandars explicitly demonstrated that, even if the 
electron has a non-zero EDM, the EDM of an atom in the non-relativistic case will be zero [4]. 
However, he went on to show that, in the relativistic treatment of the problem the atomic EDM 
does not vanish [4, 8]. Incidentally, the paramagnetic atoms, namely those with one valence elec- 
tron, are predominantly sensitive to the electron EDM [1, 2]. The electron EDM contribution 



to atomic EDM in paramagnetic atoms, is shown to be proportional to d e a 2 Z 3 [1, 2], where, 
Z is the atomic number, a is the fine structure constant and d e is the electron EDM given by, 
d e = d e (5 a where, d e is the magnitude of the electron EDM, (3 is the Dirac matrix and a is the 
spin of the electron. The atomic EDM points in the direction of total angular momentum of 
the atom. Thus, the atomic EDM due to the electron EDM would be significantly enhanced in 
heavy paramagnetic atoms. 

Indeed, we calculate a dimensionless quantity called enhancement factor (R) which is the ra- 
tio of an atomic EDM (D a ) to electron EDM (d e ) using a highly correlated all-order relativistic 
many-body theory known as the relativistic coupled-cluster (RCC) theory. In this paper, we 
present R due to the electron EDM of the ground states of Rb and Cs atoms. The current best 
limit for the electron EDM comes from an experiment on atomic thallium (Tl) [9] and prior to 
that an experiment on atomic Cs had set the limit for the electron EDM [10]. New experiments 
on the EDM of atomic Rb and Cs using the modern techniques such as laser cooling and trap- 
ping are underway in two different laboratories in the world [11, 12]. In combination with the 
results of our enhancement factors reported in this paper, they could improve the limit of the 
electron EDM by at least two orders of magnitude. 

In the following section, we have briefly described the RCC method used in calculating the 
enhancement factors of Rb and Cs atoms. In section 3, we have presented the results and com- 
pared them with the existing calculations. In section 4, we have drawn the conclusions based 
on our results and described their implications for the particle physics. 



2. RCC Method 

At first, we obtain a reference wave function (|$o)) f° r the closed-shell configuration with N — 1 
electrons, N is the total number of the electrons in the system including a valence electron (v), by 
solving the Dirac-Hartree-Fock (DF) equations with the Coulomb interaction; the Hamiltonian 
is given by, 

H = J2i Ca i ■ Pi + (& " !) m ° 2 + V nuc (n)} + zZ V c M W 

i i<j 

where, c is the speed of light in vacuum, a and (3 are the Dirac matrices, V nuc is the nuclear 
potential. 

In the framework of the RCC theory, we construct the exact wave function (|^o)) f° r the 
closed-shell configuration as, 

l*o> = e Tm |*o>, (2) 

where, are the excitation operators from the core electrons; the details of which are given in 
[13]. It is the sum of all single, double, triple and multiple excitations of occupied electrons. For 
a single valence open-shell atomic system, which is of our interest, we construct the new DF wave 
function \$> v ) by appending the valence electron v to the closed-shell DF wave function |<J>o) an d 
obtain the corresponding reference state; i.e. = ajj^o) where, a\ is the creation operator. 
The exact wave function for the corresponding valence electron system can be expressed as, 



(3) 



where, S v °^ corresponds to the excitation operator for the valence and valence-core orbitals. 
Since the systems considered in this paper have only one valence electron, the non-linear terms 
in Sy will not exist and the above wave function reduces to the form [14], 

|^) = e r(0) {l + S(°)}|^>. (4) 

It is impractical to consider all the correlated excitations in our calculation because of the 
huge requirement of the computer memory. It has been found that the RCC theory with both 
single and double excitations is quite successful in incorporating the maximum correlation ef- 
fects because of its all-order nature. Hence, we have used in this calculation, the CC method 
with only the single and the double excitations (CCSD method); i.e, = Tl^ + T2^ and 
Sv ^ = Sl v 0) + S2£ 0) . However, we have also considered the non-linear terms up to (Tl(°)) 4 . 
In our earlier works, we have given the details of the working equations to determine both the 
T( ) and S v °^ amplitudes and the ionization potential for the corresponding valence electron v 
[14]. We have also considered the leading triple excitations using CCSD(T) approximation, the 
details of which can be found in [15], to improve the accuracy of the results. 

The unperturbed RCC operator amplitudes, and S v °^ are solved in two steps: first, we 
solve for the closed-shell amplitudes, using the following equations; 



<$o|^|$o> = A£ corr (5) 
<<&S|fl?>o> = 0, (6) 

where, AE corr = Eg — E^p is the correlation energy, which is the difference between ground state 
energy and DF energy, for the closed-shell system, |3>q) represents all possible singly and doubly 

excited states with respect to the reference state, |$o) an d we define H$f = e~ T(0) H{f e T<0) 
where, H^f is the normal ordered atomic Hamiltonian given by Hq' = H$ — ($o|-^o|^o)- The 
above two equations are solved simultaneously to get T^ 1 amplitudes self-consistently. 

Then, we solve the following two equations to obtain S v °^ amplitudes and the ionization 
potential (IP) of the valence electron; 



(<S> v \{H{f} op {l + S^}\$ v ) = -AE V (7) 
{K\{WU{l + Si 0) }\*v) = -(*M 0) \*v)AE v , (8) 

where, AE V is the negative of the IP of the valence electron v and |3>*) represent all possible 
singly and doubly excited states with respect to the state | <!>„). Here, {H^} op represent all the 

open contracted operators from H^f with N electrons. The fully contracted terms are subtracted 
from the total energy E v to get AE V . 

In the absence of any external field, we have intrinsic electron EDM contribution to atomic 
EDM given by, d e Eint where, E int = -V • V nuc (n) + J2i<j Vc (nj)} is the total 

internal electric field in an atom. The energy shift due to the EDM can only be measured in 
the presence of an externally applied electric field. The applied field also induces the EDM to 
an atom which is given by, J2i{ e ?i + d e [3ai\ ■ E, where, E is the applied electric field. 

Thus, the total atomic Hamiltonian in the presence of EDM as a perturbation is given by, 



H = Hq + Hedm 



(9) 



where, H$ is the unperturbed Hamiltonian given by eqn. (1) and Hedm is the EDM perturbed 
Hamiltonian given by, 

Hedm = - d e ~^(3 cti • E int . (10) 

i 

The EDM perturbed Hamiltonian will effectively reduces to the one-body operator given by, 

HE f DM = 2icd e J2l5f3pi 2 (11) 

i 

where, f3 and 75 are the Dirac matrices, pi is the four-momentum of the electron, Ti is the 
modified Planck's constant (i.e, J^). The Hedm mixes the atomic states of opposite parities 
but with the same angular momentum. As its strength is sufficiently weak, we consider only up 
to the first-order perturbation in wave function and the modified atomic wave function for the 
valence electron state V is given by, 

|<) = |^°)) + d e |^ 1 )) (12) 

In the RCC ansatz, the cluster operators for calculating the perturbed wave function are 
given by 

T = T (0) + d e T (1) 

S v = 4°) + d e (13) 

where, and are the first order corrections to the unperturbed cluster operators 
and si°\ respectively. 

The perturbed cluster amplitudes, and are obtained by solving the following 

equations respectively in a self consistent manner; 

(^\HjT^+H e Jl M \^) = 0. 



<d>;|« - AE V ) fiW + «T« + H e Jl M ){l + Sj®}\* v ) = 

The expectation value of an atomic EDM for the state \^' v ) (with the normalized unperturbed 
wave function) is given by, 

(*' \D\W) 

a ~ <*<°>|*«°>) 

where, D = ez is the electric dipole operator. On substituting the expressions of the wave 
functions and keeping the terms only up to linear order in perturbed operators, we get the 
electron EDM enhancement factor (R= -^ £L ) as, 

{{1 + S^}DW{TV> + T^si 0) + Si 1] } + {S^ + si° V T^ + T^}W){1 + ^ 0) }1 \<f> v ) 
<d>„| e^ (0)t e^ (0) + ^ 0)t e^ (0)t e^ (0) ^ 0) |$„> 

where,m = e- Tm De TW . 



Table 1. The comparison of EDM enhancement factors of the ground state of Rb and Cs 
atoms. The references are given in the parenthesis. 



Atom 


Ground State 


EDM enhancement factor, R 


Our work 


Others 
DF + Corr. 


DF result 


DF + Corr. 


Rb 


4p b 5s (*S 1/2 ) 


19.6 


24.6 


25.67 [17] 


Cs 




94.2 


125.9 


114.2 [18] 
130.5 [19] 
120 [20] 



Table 2. The individual contributions of the important RCC terms to the EDM enhancement 
factors of Cs and Rb atoms. 



RCC term + c.c. 


Contribution (Cs) 


Contribution (Rb) 




5.33 


1.05 




138.03 


26.08 


S2( 1 ) t D(«) 


-7.15 


-1.02 


< S , l( 1 ) t D(°)5l(°) 


-2.28 


-0.39 


5l( 1 ) t £»(°)52(°) 


-6.30 


-0.95 




-0.57 


-0.07 


52(°) + L>(°)ri( 1 ) 


0.63 


0.13 


S2(°) t DS2( 1 ) 


1.97 


0.26 




0.32 


0.04 


Rest of the terms 


-1.13 


-0.14 


Normalization 


-2.91 


-0.42 


Total value 


125.94 


24.57 



3. Results & discussion 

We have computed the single particle energies and wave functions self consistently using the 
Dirac-Fock (DF) method [16] and calculated the EDM enhancement factors as explained above 
for the ground states of Rb and Cs atoms. The Table- 1 summarizes these results. The DF result 
is based on the independent particle model (IPM) approach where the electron correlation ef- 
fects are neglected. We have included all possible single, double and partial triple excitations of 
all-order using a highly correlated many-body method called the RCC theory which supplements 
the missing correlations in DF theory. 

In order to emphasize the contributions from the correlation effects to the total result, we 
have shown explicitly the individual contributions from the important RCC terms in the present 
calculation, in Table-2. It is evident from our results that, the electron correlation effects are 
quite important in the calculation of the EDM enhancement factors for the ground state of both 
Rb and Cs atoms. 

Our RCC calculation of the enhancement factors are in reasonable agreement with the pre- 



vious calculations [17, 18, 19, 20]. The calculations by Shukla et al. [17] and Das [19] are based 
on a hybrid method combining the salient features of the relativistic many-body perturbation 
theory (MBPT) and the multi-configuration Dirac-Fock (MCDF) method. Hartley et al. have 
used a variant of the relativistic MBPT [18] to perform their calculation. Their treatment of 
pair correlation effects like the former two calculations is not as comprehensive as ours. Sandars' 
calculation is based on a simple relativistic one electron theory [20]. 

As mentioned before, our calculated enhancement factors along with the results of the EDM 
experiments on atomic Rb and Cs when they come to fruition could provide the most sensitive 
limit for the electron EDM. 

There are many non-standard models (non-SM) like the different super-symmetric models, 
left-right symmetric model etc., which predict the EDM of the electron to be almost in the 
same range, where as, the SM prediction is ten to twelve orders of magnitude less which is far 
off from the reach of current experimental facilities [21]. If one can unambiguously observe the 
non-zero EDM of Rb and Cs atoms then one can combine it with our present calculation of the 
enhancement factors and choose the correct model of particle physics which predict the value of 
electron EDM in that range. It also constrains the parameter space of many models of CP vio- 
lation. Thus, the work on electron EDM has the potential to probe a new physics beyond the SM. 



4. Conclusion 

We have calculated the electron EDM enhancement factors of two paramagnetic atoms such 
as rubidium and caesium in their ground states using a highly correlated all-order relativistic 
coupled-cluster theory. Our results can be used in combination with those of the proposed high 
precision EDM experiments to obtain the stringent limits for the electron EDM, there by shed- 
ding light on leptonic CP violation which hopefully will unveil a new arena of physics beyond 
the most celebrated model of particle physics till date, the standard model. 
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